Abstract. The ISOPHOT Serendipity Survey (ISOSS) 
covered approximately 15 % of the sky at a wavelength of 
170 /im while the ISO satellite was slewing from one tar- 
get to the next. By chance ISOSS slews went over many 
solar system objects (SSOs). We identified the comets, 
asteroids and planets in the slews through a fast and ef- 
fective search procedure based on N-body ephemeris and 
flux estimates. The detections were analysed from a cal- 
ibration and scientific point of view. Through the mea- 
^vi ' surements of the well-known asteroids Ceres, Pallas, Juno 
. and Vesta and the planets Uranus and Neptune it was 
' possible to improve the photometric calibration of ISOSS 
, and to extend it to higher flux regimes. We were also able 
5h ' to establish calibration schemes for the important slew 
^ . end data. For the other asteroids we derived radiometric 
1 diameters and albedos through a recent thermophysical 
model. The scientific results are discussed in the context of 
our current knowledge of size, shape and albedos, derived 
from IRAS observations, occultation measurements and 
lightcurve inversion techniques. In all cases where IRAS 
observations were available we confirm the derived diam- 
eters and albedos. For the five asteroids without IRAS 
£SJ , detections only one was clearly detected and the radio- 
ed ' metric results agreed with sizes given by occultation and 
HST observations. Four different comets have clearly been 
detected at 170 ^m and two have marginal detections. The 
, observational results are presented to be used by ther- 
' mal comet models in the future. The nine ISOSS slews 
over Hale-Bopp revealed extended and asymmetric struc- 
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O ' tures related to the dust tail. We attribute the enhanced 



emission in post-perihelion observations to large particles 



^ ' around the nucleus. The signal patterns are indicative of 
• • , a concentration of the particles in trail direction. 

> ■ 

, Key words: Minor planets, asteroids - Comets: general 
!h ' ~ Planets and satellites: general - Infrared: Solar system 
. 5r i - Standards - Surveys 



A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

missing; you have not inserted them 



ASTRONOMY 

AND 
ASTROPHYSICS 



Solar System Objects in the ISOPHOT 170 /xm Serendipity 
Survey * 

T. G. Miiller 1 2 , S. Hotzel 3 , and M. Stickel 3 

1 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrafie, 85748 Garching, Germany 

2 ISO Data Centre, Astrophysics Division, Space Science Department of ESA, Villafranca, P.O. Box 50727, 28080 Madrid, 
Spain (until Dec. 2001) 

3 ISOPHOT Data Centre, Max-Planck-Institut fur Astronomic, Konigstuhl 17, 69117 Heidelberg, Germany 
Received / Accepted; compilation date: February 1, 2008 



1. Introduction 

The Infrared Space Observatory (ISO) (Kessler et al. 



1996) made during its lifetime between 1995 and 1998 



more than 30 000 individual observations, ranging from 
objects in our own solar system right out to the most dis- 
tant extragalactic sources. The solar system programme 
consisted of many spectroscopic and photometric stud- 
ies of comets, asteroids, planets and their satellites at 
near- and mid- infrared (near-/mid-IR) wavelengths be- 
tween 2.5 and 45 fim. At far-infrared (far-IR) wave- 
lengths, beyond 45 /im, the programmes were limited to 
spectroscopic observations of the outer planets, 3 satel- 
lites (Ganymede, Callisto, Titan), 3 comets (P/Hale- 
Bopp, P/Kopff, P/Hartley 2) and 4 asteroids (Ceres, 
Pallas, Vesta, Hygiea). Far-IR photometry on solar sys- 
tem objects was mainly done for calibration purposes 
(Uranus, Neptune and a few asteroids) and scientific stud- 
ies of extended sources (P/Hale-Bopp, P/Kopff, P/Wild 
2, P/Schwassmann-Wachmann, Chiron, Pholus). 

Additionally to the dedicated programmes on individ- 
ual sources, ISO also made parallel and serendipitous ob- 
servations of the sky. ISOCAM observed the sky in par- 
allel mode a few arc minutes away from the primary tar- 
get at wavelengths between 6 and 15 /im (Siebenmorgen 
et al. |2000| ). The ISOPHOT Serendipity Survey (Bogun 
et al. 1996) recorded the 170 /im sky brightness when the 
satellite was slewing from one target to the next. LWS per- 
formed parallel (Lhn et al. 2000) and serendipitous surveys 
(Vivares et al. 2000) in the far-IR. These complementary 
surveys contain many interesting objects, but the scien- 
tific analysis only began recently. The two LWS surveys 
and the ISOCAM parallel survey just underwent a first 
data processing and could therefore not be considered in 
the following. 
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Source extraction methods for the ISOPHOT 
Serendipity Su rvey (I SOSS) have been developed by 
Stickel et al. (1998a, 1998b) for point-sources and by 
Hotzel et al. (2000) for extended sources. First scientific 
results were published recently (Stickel et al. 



2000 



^ Toth 

et al. |2000j ; Hotzel et al. [2001D and work is ongoing to 
produce further catalogues of Serendipity Survey sources. 
To facilitate the production of source lists it is necessary 
to identify and exclude all SSOs from the survey data. 
This catalogue cleaning aspect was one motivation for the 
following analysis, but there are also the calibration and 
scientific aspects of the SSO investigations: A few well 
known asteroids and planets like Uranus and Neptune 
provide the possibility to test and extend the photometric 
calibration of ISO SS to higher brightness levels (Miiller 
& Lagerros 1998| ). For asteroids with known diameters, 
the surface regolith properties can be derived from the 
emissivity behaviour in the far-IR where the wavelength 
is comparable to the grain size dimensions. Additionally, 
reliable far-IR fluxes of asteroids allow diameter and 
albedo determinations for the less well-known targets. 
For comets the far-IR information is useful for coma and 
tail modeling (Grim et al. 2001 ). The close connection 
between large particles and far-IR thermal emission 
also allows further studies of trail formation and the 
important processes of dust supply for the interplanetary 
medium. 

In the following sections we present and discuss 
the ISOSS data with emphasis on solar system targets 
(Sect. ||). This also includes the data processing, point- 
source extraction and calibration aspects. An iterative and 
fast method to search for SSOs in large sets of slewing 
data is explained in Sect. |j. The encountered SSOs are 
then separated into 2 categories: 

1) Well known asteroids and the planets Uranus and Nep- 
tune, which were used to test and extend the photometric 
calibration of ISOSS (Sect. @). 

2) Asteroids and comets, for which the far-IR fluxes were 
used a) to derive diameters and albedos (asteroids) or b) 
to give a qualitative and quantitative description of the ob- 
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Fig. 1. IRAS 100 /im map with the slew paths of the 4 
C200 pixels overplotted. The symbols used for the indi- 
vidual detector pixels are the same as in Fig. ^. Here, 
their colour coding indicates qualitatively the measured 
intensities. The intensity peaks in Pixel 1 and 4 have no 
correspondence in the IRAS map. The star symbol marks 
the N-body ISO-centric position of Ceres at the time of 
the slew. Measurements within 5' from the detector cen- 
tre are marked here and in Fig. || as black crosses. The 
detector and pixel apertures as well as the scan direction 
are indicated. The slew TDT number is 09380600 (see also 
Table B>. 



servational results for future modeling (comets) (Sect. ||). 
In Sect. ^ we summarize the results and give a short out- 
look to future projects. 



2. Observations and Data Reduction 

ISOPHOT Serendipity Survey (ISOSS) measur ement s 
were obtained with the C200 detector (Lemke et al. 1996| ), 
a 2x2 pixel array of stressed Ge:Ga with a pixel size of 
89.4". A broadband filter (C_160) with a reference wave- 
length of 170 (an and a width of 89 was used. The high- 
est slewing speed of the satellite was 8'/sec. During each 
1/8 sec integration time 4 detector readouts were taken, 
i.e. the maximum read out distance on the sky was 15" 
yielding one brightness value per arcminute (see Figs. [I] 
and |). During the ISO lifetime, about 550 hours of ISOSS 
measurements have been gathered, resulting in a sky cov- 
erage of approximately 15%. 
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Fig. 2. The calibrated pixel intensities as a function of 
slew length. Intensities of Pixel 2-4 are shifted downwards 
in steps of 20 MJysr -1 . The slew section shown corre- 
sponds to Fig. |l| Black crosses mark measurements with 
detector positions closer than 5' to the calculated asteroid 
position. Pixel 1&4 cross Ceres almost centrally. 



2.1. Data Analysis 

A standard data processing was applied using the 
ISOPHOT Interactive Analysis PIAQ (Gabriel et al. [1997| ) 
Version 7.2 software package. The detailed processing 
steps are given in Stickel et al. ( |200Cl) ) . Special care had 
to be taken to correct gyro drifts between sequent guide 
star acquisitions of the star tracker. For point-sources, 
the deglitched and background subtracted signals of the 
4 pixels were phase-shifted according to the position an- 
gle of the detector and co-added. Source candidates were 
searched for in this co-added stream by setting a cut of 3 a 
of the local noise. Then, the source position perpendicu- 
lar to the slew was determined from a comparison between 
signal ratios with a gaussian source model. The flux was 
afterwards derived from 2-D gaussian fitting with fixed 
offset position. 

In case of long slews, the surface brightness were de- 
rived from a measurement of the on-board Fine Calibra- 
tion Source (FCS) preceding the slew observation. For 
short slews the default C200 calibration was used. To 
tie point-source fluxes derived from ISOSS to an abso- 
lute photometric level, dedicated photometric calibration 
measurements of 12 sources, repeatedly crossed with vary- 
ing impact parameters were compared with raster maps 



on the same sources (Stickel et al. 1998a). The compar 



1 The ISOPHOT data presented in this paper were reduced 
using PIA, which is a joint development by the ESA As- 
trophysics Division and the ISOPHOT Consortium with the 
collaboration of the Infrared Processing and Analysis Cen- 
ter (IPAC). Contributing ISOPHOT Consortium institutes are 
DIAS, RAL, AIP, MPIK, and MPIA. 
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ison between slew and mapping fluxes showed that for 
brighter sources the slewing observations miss some signal, 
probably due to transient effects in the detector output 



(Acosta-Pulido et al. 2000) in combination with detector 
non-linearities. For sources brighter than 30 Jy Stickel et 
al. ( 2000D found signal losses of 50%, although the true 
losses were not well established due to a lack of reliable 
sources. For fainter sources (<10Jy) the flux loss in the 
slews is only 10-20%. 

2.2. Source Extraction Methods 

The SSOs were encountered at different slew speeds, which 
can be characterized by 'fast': above 3'/sec; 'moderate': 
1.5'/sec < speed < 3'/sec; 'slow': below 1.5'/sec; 'stop': 
at slewends, like a staring observation. But aspects like 
the background level, the detector history and impact pa- 
rameters also play a crucial role in source extraction and 
flux calibration methods. 

2.2.1. Method 1 



Automatic point source extractor (Stickel et al. 2000 ) 
for all slewing speeds above 1.5'/sec and non-saturated 
crossings. All source candidates were cross correlated with 
the list of SSO candidates (see Fig. ||) and the associa- 
tions found carefully examined. Flux loss corrections (see 
Sect. 2.1) have to be applied. 



2.2.2. Method 2 

For all slewing speeds, but using only the pixel with the 
highest signal and converting it to flux density as if the 
source was centred. This method leads to upper and lower 
flux limits only. The lower limits, designated by '>' or 
';§>', are connected to clear detections. The quality of the 
lower estimate depends on the impact parameter. Useful 
upper limits have only been given for direct hits where 
no detection signal was seen. The upper limit then corre- 
sponds to the 3cr- valu e of the background noise. Flux loss 
corrections (Sect. ^l|) have to be applied as for Method 1. 
Note: It is assumed that 64% of the flux density of a pixel 
centred point source falls onto this pixel. 



2.2.3. Method 3 

This method has been used at slewends, if the source was 
inside the detector aperture: 

a) Using the signals of all 4 pixels at the very end of the 
slew and converting them to flux densities assuming 
the source is centred on the detector. Note: Only 53% 
of the flux of an array centred point source is detected, 
of which 21% fall on pixel 1, 24% o n pixe l 2, 32% on 
pixel 3 and 23% on pixel 4 (Laureijs 1999 ). 



b) Like case 3a, but source centred on one pixel. Note: In 
total 74.3% of the source flux are seen by the 4 pixels: 
64% in the source pixel, 2x4.2% in the two adjacent 
pixels and 1.9% in the diagonal pixel. 

The statistical errors are computed from the weighted re- 
sults of the 4 pixels. 

Not all of the ISO scientific targets are to be found in 
the end-of-slew data: Firstly, the 4 ISO instruments view 
separate areas of the sky. Slew end position (ISOPHOT) 
and target position (other instrument) can therefore differ 
by up to 20'. Secondly, many observing modes, especially 
for ISOPHOT, started off-target for mapping purposes or 
to avoid strong detector transients. 

3. Solar System Object Identification 

The identification and separation of moving solar sys- 
tem targets from the Serendipity slews is difficult: The 
Serendipity slew data consist of very narrow stripes across 
the sky, lacking, to first approximation, any redundancy. 
Additionally, the colour information is missing and in the 
far-IR region the cirrus confusion is a serious problem. 
Therefore, the solar system object identification was done 
on basis of accurate ephemeris calculations and model flux 
estimates (see Miiller [2001 ). 

The Serendipity slew data consist of 11 847 slews with 
a total length of 141 411°. The slews were cut in 4 232 525 
individual pointings of approximately 2' length. Each of 
these pointings had to be checked against SSOs. On 20th 
of March 2000, the Minor Planet Center archive con- 
sisted of 68840 asteroids (14308 numbered, 24598 unnum- 
bered with multiple-opposition orbits and 29934 unnum- 
bered with single-opposition orbits) and 237 comets. Ad- 
ditionally, the outer planets and their satellites had to 
be included, leading to a total of approximately 3 • 10 11 
ephemeris calculations. It was therefore necessary to pre- 
select the number of SSOs considerably and to invent fast 
search procedures. 

3.1. SSO Preselection 

To facilitate and speed up the search process, only SSOs 
which at maximum are brighter than the sensitivity limit 
of 1 Jy at 170 /Ltm have been considered: The outer plan- 
ets (Mars, Jupiter, Saturn, Uranus, Neptune and Pluto) 
and their satellites were included. The inner planets were 
not visible for ISO. Due to the difficulties of predict- 
ing the brightness of active comets, no initial flux pres- 
election was done for the 237 comets. The unnumbered 
single-apparition asteroids have not been considered be- 
cause of possible large ephemeris uncertainties and gen- 
erally too low brightness at 170 ^m (except for a few 
Near-Earth asteroids which can reach this flux limit at ex- 
tremely close encounters). The preselection of numbered 
and unnumbered multi-apparition asteroids was based on 
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Fig. 3. Extraction procedure for ISOSS slew data to find 
SSO candidates, based only on pointings and timings of 
ISO slews and expected 170 /im fluxes. The flux preselec- 
tion is described in the text. 



conservative flux calculations, using a simplified Stan- 
dard Thermal Model (STM, Lebofsky et al. |1986| , |1989[ ) 
and assuming a non-rotating spherical object. The fol- 
lowing conservative input values were used: p v = 0.02, 
G = —0.12, e = 1.0, f] = 0.7. The input diameter was 



calculated using p v and the absolute magnitude H with: 
log p v = 6.259-2 log D eff -0AH. In Flux Preselection I, 
with a cut limit of 1.0 Jy, the object was assumed to be 
located at perihelion during opposition. This reduced the 
number of asteroids by 98% from 38906 to 596. In Flux 
preselection II, with a cut limit of 0.5 Jy, the real calcu- 
lated distances (r, A) were used. 

For the comets we determined 36 "hits" where the ob- 
ject was within 5' from the slew centre and with the comet 
being within 3AU (for Hale-Bopp 5AU) from the Sun. 
A simple flux estimate lowered the number to 16 possi- 
ble candidates. The 170 fim flux estimate was based on 
an assumed dust albedo A = 0.10 and a temperature of 
dust particles of To = 330 K at the heliocentric distance 
of r = 1 AU through the following formula: 

F v = B(X,T) ■ f ■ ^ (1) 



with the temperature T = ^/T 4 • (1 - A)/r 2 and A the 
geocentric distance to the comet. We assumed a dense 
central coma of 10 000 km radius with a filling factor of 
/ = 10~ 4 and 1/R brightness profile out to a distance 
of 50 000 km. The corresponding model predictions fitted 



nicely the published results by Campins et al. 1990 for 



comet P/Tempel 2 and by Hanner et al. 1994 for comet 
Mueller 1993a. This model was also used for initial flux 
estimates for the preparation of ISO comet observations 
(Grim, private communication). 



3.2. Search Radius 

The search radius for each pointing had to be much bigger 
than the real 3' x 3' field of view of the detector for several 
reasons: 1) slew position uncertainties (up to 2', Stickel 
et al. 2000| ); 2) uncertainties of the 2-body unperturbed 
ephemeris, based on 200-day epoch orbital elements (up 
to a few arcmin); 3) the ISO parallax (up to 3' for close 
encounters at 0.5 AU). In the first iteration, the search 
radii were set to 8.1' for asteroids (3' for ephemeris un- 
certainties, 3' for ISO parallax, 2.1' for the centre-corner 
distance of the C200 array), to 30' for comets (to account 
for extended structures) and to 2° for the bright plan- 
ets (to account for possible straylight influences). In the 
second iteration, after the ephemeris recalculation with 
an N-body programme and after parallax corrections, the 
search radius was uniformly set to 5'. Additionally, all 
identified slews from the first iteration were marked, be- 
cause of possible influences from bright SSOs. 



3.3. Search Procedure 

Figure ^| summarizes the procedure in detail, giving also 
the input and output number of targets. With this pro- 
cedure it was possible to reduce the initially estimated 
10 11 ephemeris calculations to 10 9 2-body and 10 4 N- 
body calculations. The final potential hits of 56 asteroids, 
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16 comets and 22 planets fulfilled the flux requirements 
at the actual time of the observation and were located 
within 5' of the slew. Note that we count each encounter 
of a slew with an SSO as "hit" . The actual numbers of dif- 
ferent objects in this list are 21 asteroids, 7 comets and 2 
planets. These results, based on pure pointing and timing 
information, are strongly influenced and biased by the un- 
equal distribution of the slews in the sky, satellite visibility 
constraints and the ISO observing programme itself. The 
relative large number of comets is due to the weak flux 
limit and it was clear that not all of them would be bright 
enough to be detected. 

4. Calibration Results 

ISOSS observations of Uranus, Neptune and well known 
bright asteroids enabled us to improve and extend the ex- 
isting ISOSS calibration (Method 1). They also allowed us 
to estabish the calibration of new source and flux extrac- 
tion methods, namely Methods 2, 3a and 3b (see Sect. ||). 

The Uranus and Neptune models are based on Grif- 
fin & Orton ( 1993| ) and Orton & Burgdorf (priv. comm.), 
respectively. For Ceres, Pallas, Ju no and Vest a athermo- 
physical model (TPM) (Lagerros |l996j ; |l997| ; |l998| ) was 
used to predict their brightnesses at the times of the obser- 
vations. The TPM and its in put parameters are des cribe d 
in Miiller & Lagerros ( |l998| ) and in Miiller et al. ( |l999| ). 
The quality and final accuracy of TPM predictions are 
discussed in Miiller & Lagerros ( 2002a ). The general as- 
pects of asteroids as calibration standa rds for IR projects 
are summarized in Miiller & Lagerros ( 2002b ). 

Photometric measurements of different astronomical 
sources can be compared on bases of colour corrected 
monochromatic fluxes at a certain wavelength or on basis 
of band pass fluxes. In this calibration section all model 
fluxes have been modified by an "inverse colour correc- 
tion" in a way that they correspond to ISOSS band pass 
measurements of a constant energy spectrum (yF v = 
const.). This implied inverse colour correction terms of 
1.09 -1 for Uranus and Neptune (both have temperatures 
at around 60 K at 170 /im) and 1.17 - 1 for the bright main- 
belt asteroids (assumed far-IR temperature of 180-200 K), 
see also the colour correction tables in "The ISO Hand- 



book, Volume V", Laureijs et al. (2000) 



4.1. Method 1 

ISOSS crossings over planets and asteroids, which were de- 
tected by the Automatic Point Source Extractor (Stickel 
et al. 2000), are listed in Table [l], where the columns are: 
(1) TDT number of the slew, (2) date and Universal Time 
at the moment of the SSO observation, (3) name of the so- 
lar system object, (4) observed flux density, (5) predicted 
flux density, (6) ratio between observed and modeled flux 
density (see also Fig. |). The ISOSS results are the FCS 
calibrated band fluxes. The model predictions were modi- 
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Fig. 4. The ratio of Serendipity slew flux densities and 
model predictions for reliable Uranus, Neptune, Ceres, 
Pallas, Juno and Vesta observations. For bright sources, 
the Serendipity slews miss some flux. 



ficd by an inverse colour correction to make them compa- 
rable with the ISOSS measurements (see above). All list 
entries of Uranus, Neptune, Ceres, Pallas and Vesta give 
a ratio between observed and model flux of (0.58 ± 0.05), 
for fluxes larger than about 25 Jy. At fluxes below 25 Jy 
(only 2 cases) the ISOSS to model ratios are close to 1.0. 
This is in ex cellent agreement with the results of Stickel 
et al. J2000[ ). They showed that ISOSS slew fluxes of 12 
selected galaxies were systematically lower than fluxes de- 
rived from dedicated maps. To bring the fluxes from map- 
ping and slewing into agreement ISOSS fluxes larger than 

30 Jy were corrected with an estimated constant scal- 
ing factor of 2, while lower fluxes were scaled with a flux 
dependent correction function. Table |l| represents there- 
fore the first direct flux calibration of the PHT Serendipity 
Mode as compared to the previously used indirect method 
of flux ratios between PHT22 raster maps and slew results. 

Figure |^ shows the ratios between the flux densities 
derived from ISOSS and the 170 fim model predictions. 
The stars represent the results from dedicated calibration 
measurements (Stickel et al. 2000D , the filled circles are 
values from Table [l] Uranus, Neptune, Ceres, Pallas, Juno 
and Vesta, serendipitously seen by ISOSS, provide now a 
reliable calibration at higher flux densities. 



4.2. Method 2 

Table ^ summarizes the values which were derived from 
the solar system far-IR standards for slow slewing speeds, 
saturated measurements and sources outside the slews. 
These measurements were rejected by the source extrac- 
tion procedures of Method 1. The table columns are: (1-6) 
same as in Table [j], (7) slew speed category at the moment 
of the SSO observation, (8) additional remarks. 
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TDT Date/Time SSO F 0bs F Model F 0bs /F Model 

No. (Jy) (Jy) 



(1) 








(2) 


(3) 


(4) 


(5) 


(6) 


07881200 


03-FEB-96 


09 


:46: 


:42 


(4) Vesta 


28.9 


39.6 





.73 


10180400 


26-FEB-96 


06 


:20: 


:00 


(4) Vesta 


30.7 


54.2 





.57 


14080700 


05-APR-96 


15 


:39: 


:10 


Neptune 


145.1 


271.3 





.53 


23080100 


03-JUL-96 


20 


:15: 


:13 


(2) Pallas 


15.2 


27.5 





.55 


32181100 


03-OCT-96 


04 


:38: 


:13 


Neptune 


153.2 


279.8 





.55 


42283300 


ll-JAN-97 


22 


:47: 


1 1 


(3) Juno 


12.6 


12.0 


1 


.05 


34480700 


26-OCT-96 


00 


:25: 


:16 


Neptune 


159.1 


272.7 





.58 


69880600 


13-OCT-97 


21 


:27: 


:19 


Neptune 


160.7 


277.6 





.58 


70681100 


22-OCT-97 


02 


:44: 


:35 


Neptune 


158.0 


274.9 





.57 


71381000 


29-OCT-97 


05 


:06: 


:45 


Neptune 


168.3 


272.7 





.62 


71980500 


03-NOV-97 


22 


:46: 


:18 


Neptune 


149.9 


271.3 





.55 


72081500 


05-NOV-97 


01 


:19: 


:05 


Uranus 


395.7 


672.8 





.59 


72081600 


05-NOV-97 


01 


:57: 


:38 


Neptune 


156.1 


270.5 





.58 


76280400 


16-DEC-97 


13 


:22: 


:05 


(1) Ceres 


31.5 


52.4 





.60 


79781500 


21-JAN-98 


00 


:30: 


:12 


(4) Vesta 


24.2 


23.9 


1 


.01 



Table 1. Results for Method 1. The model fluxes are multiplied by 1.09 (planets) and by 1.17 (asteroids) to account 
for the spectral shape differences between vF v — const, (assumed spectrum in the ISO calibration) and the real object 
spectrum. Column (4) contains the FCS calibrated fluxes. 



TDT 


Date/Time 


SSO 


^Obs 


^Model 


fbbs / -FModel 


Slew speed 


Remarks 


No. 






(Jy) 


(Jy) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


09380600 


18-FEB-96 15:10:15 


(1) Ceres 


>46 


73.8 


>0.62 


slow 


ok 


29280600 


04-SEP-96 00:31:59 


(1) Ceres 


>54 


67.2 


>0.80 


slow 


very high bgd. 


32880600 


09-OCT-96 23:11:09 


Neptune 


>155 


277.8 


>0.56 


slow 


ok 


36381700 


14-NOV-96 04:34:23 


Neptune 


>50 


267.2 


>0.19 


moderate 


outside 


54480800 


13-MAY-97 12:58:55 


Uranus 


>236 


700.1 


>0.34 


moderate 


saturated 


55280300 


21-MAY-97 06:17:09 


Uranus 


>70 


709.5 


»0.10 


stop 


outside 


69880200 


13-OCT-97 17:39:49 


Uranus 


>265 


700.1 


>0.38 


moderate 


saturated 


69880500 


13-OCT-97 20:48:47 


Uranus 


>245 


700.1 


>0.35 


moderate 


saturated 


71480300 


29-OCT-97 23:34:32 


Uranus 


>202 


680.8 


>0.30 


moderate 


saturated 


87481000 


07-APR-98 14:36:41 


Uranus 


>232 


652.0 


>0.36 


moderate 


saturated 



Table 2. Results for Method 2. The model fluxes are multiplied by 1.09 (planets) and by 1.17 (asteroids) to account 
for the spectral shape differences between vF v = const, (assumed spectrum in the ISO calibration) and the real object 
spectrum. The values in Col. (4) are already corrected for the individual pixel point-spread function (F ps f = 0.64). 



The results from Method 2 show that also difficult slew 
data with either saturated pixels, objects slightly outside 
the array or slow speeds can be used to derive useful lower 
limits for interesting sources. As the satellite still moves 
the flux loss corrections from Method 1 have to be applied 
to get the best lower limits. In fact, for the 2 unprob- 
lematic hits (TDT 9380600 and 32880600) with neither 
saturated signals nor large impact parameters, the flux 
loss correction brings the ISOSS fluxes within 10% of the 
model predictions. 

4.3. Method 3 

At the slewend, when the satellite does not move anymore, 
the ISOSS data can in principle be treated as normal C200 



photometric data. Two ideal cases - source centred on 
the array (Method 3a) and source centred on one pixel 
(Method 3b) - can be distinguished. The results on the 
bright sources for both methods are summarized in Ta- 
ble ^, where the columns are the same as in Table |. The 
uncertainties in the table, given in brackets, are statisti- 
cal errors of weighted results from all 4 pixels. The results 
of Method 3 are compared with the model predictions in 
Fig. | 

The 5 Neptune measurements (Method 3a) agree 
nicely with the model predictions (Observation/Model: 
0.96±0.09). For the fainter asteroids the Method 3a over- 
estimates the flux systematically by 10-50%, depending on 
the brightness level (see Fig. |j) . The discrepancy between 
bright and faint sources is probably due to detector nonlin- 
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TDT Date/Time SSO F 0bs ^Modd F 0hs /F Modcl 

No. (Jy) (Jy) 



(1) 








(2) 


(3) 


(4) 


(5) 


(6) 


09380500 


18-FEB-96 


07 


11 


04 


(1) Ceres 


70.1(6.9) 


74.4 





94 


15480200 


19-APR-96 


03 


42 


10 


Neptune 


269(11.1) 


275.3 





98 


23781000 


ll-JUL-96 


04 


13 


16 


(3) Juno 


12.9(2.6) 


9.5 


1 


35 


25180400 


24-JUL-96 


23 


47 


09 


(2) Pallas 


26.4(2.0) 


22.3 


1 


18 


26580800 


08-AUG-96 


02 


50 


48 


(2) Pallas 


25.1(2.4) 


19.5 


1 


29 


27580200 


18-AUG-96 


02 


17 


46 


(1) Ceres 


85.2(10.2) 


79.4 


1 


07 


32880500 


09-OCT-96 


21 


54 


26 


Neptune 


236(18.3) 


277.8 


o 


85 


35680200 


06-NOV-96 


20 


07 


37 


Neptune 


267(16.1) 


269.3 





99 


38781200 


07-DEC-96 


23 


52 


44 


(3) Juno 


22.4(5.4) 


15.6 


1 


44 


41980900 


08-JAN-97 


18 


53 


29 


(3) Juno 


16.3(0.5) 


12.5 


1 


30 


51080600 


09-APR-97 


10 


54 


42 


(2) Pallas 


14.1(2.1) 


10.8 


1 


30 


51080800 


09-APR-97 


15 


19 


52 


(2) Pallas 


16.1(1.3) 


10.5 


1 


53 


51380100 


12-APR-97 


04 


33 


19 


(2) Pallas 


16.7(3.9) 


11.5 


1 


45 


53980100 


08-MAY-97 


03 


39 


36 


Neptune 


282(70.4) 


280.8 


1 


00 


53980300 


08-MAY-97 


11 


11 


10 


Neptune 


278(31.1) 


280.8 





99 


54581400 


14-MAY-97 


10 


49 


26 


(1) Ceres 


67.2(4.3) 


55.5 


1 


21 


57581500 


13-JUN-97 


13 


53 


04 


(4) Vesta 


32.5(2.9) 


24.6 


1 


32 


74881000 


03-DEC-97 


02 


21 


54 


(1) Ceres 


67.1(3.5) 


58.3 


1 


15 


53880300 


07-MAY-97 


07 


54 


19 


(1) Ceres 


43.8(4.9) 


52.5 





83 


61580800 


23-JUL-97 


02 


07 


07 


(4) Vesta 


31.1(1.2) 


34.4 





91 



Table 3. Results for Method 3 (upper part: 3a, lower part 3b). The model fluxes are multiplied by 1.09 (planets) and 
by 1.17 (asteroids) to account for the spectral shape differences between vF v — const, (assumed spectrum in the ISO 
calibration) and the real object spectrum. The uncertainties in the table, given in brackets, are statistical errors of 
weighted results from all 4 pixels. 



parison of Fig. || with Fig. ^ supports this explanation, as 
both diagrams show a decrease in the detector signals for 
bright sources. The fast slewing on the other hand, which 
affects Method 1 but not Method 3, could be responsible 
for the generally too low ISOSS fluxes in Fig. |^. 

Both options of Methods 3 open a powerful new possi- 
bility to evaluate the 170/im fluxes of many scientific ISO 
targets, which are quite often covered in the end of slews 
before the intended science programme started. 

4-4- Pointing Comparison 

The N-body ephemeris calculations for our SSOs included 
a transformation from geocentric to ISOcentric frame. 
The maximal geo-/ISOcentric parallax corrections were: 
737.7" for the Apollo asteroid (7822) 1991 CS, 336.6" 
comet P/Encke and 61.2" for Mars. The final accuracy 
of the ISOcentric SSO ephemeris has been estimated to 
about 1-2". 

The ISOSS signal pattern, i.e. the relative signals of 
the 4 pixels, is a very sensitive indicator of the exact po- 
sition of the source within the detector array. All close 
encounters have been checked by eye for discrepancies be- 
tween predicted slew offsets and the signal patterns. No 
disagreement was found, which implies that the predicted 
SSO positions and the slew positions agree with each other 
within 30", corresponding to 1/3 pixel width. In case of 




^Model / JY 

Fig. 5. The ratio of Serendipity flux densities from 
Method 3a and model predictions for Neptune, Ceres, Pal- 
las, Juno and Vesta observations. Error bars are statistical 
errors from the individual pixel results. Circles encompass 
data points from slews which had to be calibrated with 
the default calibration; in these cases, the true uncertain- 
ties exceed the given statistical errors. A flux dependency 
similar as in Fig. ^ (Method 1) can be seen. 



earities, which are not corrected in the OLP 7 Serendipity 
Mode data, and which could be responsible for the flux 
dependency of the scaling factor (see Sect. 4.1). A com- 
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non-detections, the SSOs were either too faint, or they 
were actually just outside the slew. This high pointing 
accuracy allowed us to give upper limits (depending on 
the background) in cases when the source was crossed by 
the slew but no signal was detected (see also Sect. 2.2.2Q . 
In slew direction the position accuracy is better than 1', 
limited by fast slewing in combination with the detector 
read-out frequency. 

5. Scientific Results 

The procedure from Sect. || resulted in a list of potential 
SSO candidates in the ISOSS. Mainly for the following 
reasons, not all of them were visible in the slew data: 

1. ) The structured and bright cirrus background caused 

source confusion and limited the point source extrac- 
tion It also affected strongly the analysis of extended 
structures, like from a cometary coma. 

2. ) The sensitivity limit of approximately Uy at 170 /mi 

allowed only the detection of bright asteroids and 
comets, which are already well known through other 
observing programmes and techniques (IRAS, occulta- 
tion measurements, radar, ...) and through dedicated 
ISO measurements. 

3. ) The source extraction from a 4-pixel camera is diffi- 

cult due to the high slewing speed and the variety of 
impact parameters. The resulting fluxes or flux limits 
had usually larger error bars than comparable pointed 
observations, where sources were usually either centred 
in the C200 array or on a single pixel. 

4. ) For some faint sources the allowed maximal offset of 5' 

was too large to produce a noticable signal increase. 

In some cases the object was visible, but a reliable flux 
determination from the ISOSS was not possible. In these 
cases upper or lower limits are given. 

Tables ||, || and || include all SSO predictions which 
are within 5' of the slew centre, fulfill the flux require- 
ments and have not been used in Sect. ||. For the sci- 
entific comparison between ISOSS fluxes and model pre- 
dictions we did the following calibration steps: We deter- 
mined the ISOSS calibrated inband fluxes through the dif- 
ferent methods and corrected them by an estimated factor 
based on the slopes visible in Figs. [| and 0. As a last step 
we applied the colour correction to obtain monochromatic 
flux densities at 170 /im (ISOSS values in Tables |, [§). 

5.1. Planets 

The inner planets were not visible for ISO due to pointing 
constraints. Mars, Jupiter and Saturn exceeded the satu- 
ration limits, Pluto was below the 1.0 Jy limit. Therefore, 

2 Even bright sources like Ceres are sometimes difficult to 
analyse if they are in regions of high background like at \ ec i. = 
90° and A ec ;. = 270°, where the ecliptic crosses the galactic 
plane. 



only Uranus and Neptune were seen, but already used to 



extend the ISOSS calibration (Sects, fll], and Q. But 
the bright planets had to be included in the search pro- 
gramme with the objective to identify close-by slews. For 
the extremely IR bright planets, diffraction effects of the 
optical system in combination with certain satellite-planet 
constellations produced bright spots, spikes and ring like 
structures around the planets, which are visible in the slew 
data. Mars, Jupiter and Saturn, with 170 /im brightnesses 
between 10 000 and 400 000 Jy, influenced slews up to 1° 
distance, Uranus and Neptune (> 200 Jy) up to 10' dis- 
tance. These slews have been identified (with the above 
mentioned SSO extraction method in combination with a 
large search radius) and the planet influence can now be 
taken into account for further scientific catalogues based 
on ISOSS. Some planetary satellites are bright enough to 
be visible in principle. However, close to Jupiter (the max- 
imal distance for the Galilean satellites is about 11') and 
Saturn (the maximal angular distance for the 8 largest 
satellites is about 10') no 170 /zm fluxes can be derived, 
due to the strong planet influences. 

The measured Uranus and Neptune flux values from 
Tables @, § and | can also be used as input to future mod- 
els of planetary atmospheres. Current models are based on 
Voyager IRIS data from 25 to 50 /im a nd su b-millimetre 
data beyond 350 /mi (Griffin & Orton |l993|) with an in- 
terpolation in between. The ISOSS observations make this 
wavelength range directly accessible for model tests in the 
far-IR. 

5.2. Asteroids 

After establishing new methods for the flux calibration, 
based now additionally on measurements of Uranus, Nep- 
tune, Ceres, Pallas, Juno and Vesta, the monochromatic 
flux densities at 170 /xm of the remaining asteroids were 
derived. 23 out of the 56 asteroid "hits" were already used 
in this calibration context in the previous Sect. | The re- 
maining 33 hits can be split in 3 groups: 

5.2.1. IRAS and Poor ISOSS Detection 

18 asteroid predictions have reliable IRAS detections, but 
only low quality ISOSS detections. The reasons for the 
poor ISOSS fluxes are manifold: slew offsets, bright back- 
grounds, technical proble ms, l ow fluxes, etc. For these 18 
asteroids Tedesco et al. ( 1992 ) calculated already diame- 
ters and albedos, based on IRAS observations. Upper flux 
limits from ISOSS would therefore not give any new in- 
formation. 



5.2.2. IRAS and Good ISOSS Detection 

7 asteroids (9 hits) have reliable IRAS detections and also 
good quality ISOSS detections (see Table [|). For these 
asteroids we could derive successfully fluxes and upper 
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TDT 


Date/Time 


SSO 




-FModel 


Method 


Slew 


Remarks 




No. 






(Jy) 


(Jy) 




speed 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




11080300 


06-MAR-96 11:42:00 


(5) Astraea 


<4 


1.3 


2 


fast 


upper limit, cirrus 


bgd. 


12780300 


23-MAR-96 02:50:17 


(5) Astraea 


<4 


1.6 


2 


fast 


upper limit, cirrus 


bgd. 


85480200 


18-MAR-98 13:48:12 


(7) Iris 


<3 


2.1 


2 


fast 


in cirrus knot 




25381400 


27-JUL-96 03:22:14 


(15) Eunomia 


2.9 


2.5 


1 


fast 


ok 




85480300 


18-MAR-98 17:22:43 


(89) Julia 


<3 


1.7 


2 


fast 


upper limit 




18780100 


22-MAY-96 05:40:00 


(344) Desiderata 


4.5 


5.8 


1 


fast 


ok 




18280800 


17-MAY-96 03:50:00 


(532) Herculina 


5.3(1.0) 


5.5 


3 a 


stop 


ok 




21780900 


21-JUN-96 06:01:48 


(532) Herculina 


3.4(0.3) 


3.8 


3 a 


stop 


ok 




83380500 


25-FEB-98 11:36:26 


(1036) Ganymed 


<4 


0.1 


2 


fast 


upper limit 





Table 4. ISOSS observational results for asteroids with IRAS detections. Note, that the ISOSS results were first flux 
corrected according to Sect. || and then colour corrected. The uncertainties in the table, given in brackets for method 
3a, are statistical errors of weighted results from all 4 pixels. 



limits with our newly established calibration, based on 
clear detections. The IRAS diameter and albedo values in 
the following are all taken from the Minor Planet Survey 
(MPS, Tedesco et al. |1992[). 



Astraea: The main-belt asteroid Astraea has been well 
observed by different techniques, including IRAS, radar 
and lightcurve observation. The combination of all mea- 
surements led to the description of the object as a rotat- 
ing e llipsoid with a well determined spin vector (Erikson 
2000|) and ax is dim ensions of 143(±12%) x 115 x 100 km 
(Magri et al. |l999|) . Using the TPM with default th erma l 
parameters for main-belt asteroids (Miiller et al. 1999 ) 
together with the shape, size and spin vector information 
gave fluxes of 1.3 ± 0.3 Jy and 1.6 ± 0.4 Jy at the 2 ISOSS 
epochs (see Table ||). The measured ISOSS upper limits 
are in agreement with the TPM predictions. 



Iris: Like for Astraea, a shape model has been established 
for Iris based on a combination of radio metric, lightcurve 
and occultation data (Magri et al. 1999 ). The correspond- 
ing TPM prediction gives 2.1 ± 0.2 Jy at the ISOSS epoch, 
with an additional lightcurve variation of about 25 % (min 
to max). The measured lower flux limit is in agreement 
with the calculations. 



Eunomia: Eunomia was one of the best observed aster- 
oids by IRAS: 7 epochs distributed over almost one month, 
each time observed with high S/N in all 4 bands. The 
MPS diameter is 255.3 ± 15 km and the albedo 0.21 ± 
0.03. Two single chord occultation measurements led to 
diamete rs of > 309 ± 5 km (Overbeek |l98|) and >232km 
(Stamm |l99l] ). The TPM prediction (based on MPS diam- 
eter and albedo and on shape and spin vector by Erikson 
2000| ) gives 2.5 ± 0.5 Jy at the ISOSS epoch, with the main 
error contribution coming from the large lightcurve varia- 



tion. The measured ISOSS flux of 2.9 Jy agrees within the 
errorbars. 



Julia: IRAS observed this asteroid 4 times within 2 
weeks, each time with high S/N in all 4 bands. The MPS 
diameter is 151 ± 3.1km and the albedo 0.18 ± 0.01. 
No shape or spin vector is available, but the possible 
lightcurve amplit udes r ange between 0.10 and 0.25 mag 
|1989|) . The TPM prediction (based on 



(Lagerkvist et al. 
MPS diameter and albedo together with a spherical shape) 
gives 1.7 ± 0.1 Jy at the ISOSS epoch, with an additional 
maximal lightcurve variation of approximately 25 % (min 
to max). The measured upper limit agrees with this pre- 
diction. 



Desiderata: This asteroid was observed by IRAS ex- 
tensively at 9 epochs during a period of 2 months with 
high S/N in all bands. The MPS diameter is given with 
132.3 ± 5.5 km and the albedo 0.06 ± 0.01. No shape and 
spin vector exists currently for Desiderata, but a 0.17 mag 
lightcurve amplitude has been stated by Lagerkvist et 
al. (p)89|). The TPM prediction (based on MPS diame- 



ter and albedo together with a spherical shape) gives 5.8 
± 0.5 Jy at the ISOSS epoch, with an additional maxi- 
mal lightcurve variation of approximately 17% (min to 
max). Assuming an ISOSS observation at lightcurve min- 
imum and a diameter at the lower end of MPS diameter 
range would result in a TPM flux which is only a few per- 
cent above the measured ISOSS value, but well within the 
ISOSS measurement error bars. 

Herculina: 7 IRAS observations with high S/N in either 
3 or 4 bands have been obtained between March and Oc- 
tober 1983. The The MPS diameter is given with 222.2 
± 7.6 km and the albedo 0.17 ± 0.01. The occultation 
diameter of 217 ± 15 km is based on several chords in 
combination with information on the pole orientation and 
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a lightcurve fit (Bowell et al. 1978). MPS and occulta- 
tion diameters agree nicely. The complete shape and spin 
vector solutions derived from lightcurve observations are 
given in Erikson ( |2000 ). Using the full information for Her- 
culina (see also Miiller & Lagerros 1998 for details) led to 
170 ^m fluxes of 5.5 ± 0.4 Jy and 3.8 ± 0.3 Jy. These values 
have been calculated using the exact lightcurve phase and 
amplitude at the time of the ISOSS observations. The al- 
most perfect agreement between predicted and measured 
fluxes (see Table confirms in an independent way the 
reliable calibration of this new flux extraction method for 
the ISOSS. 



Ganymed: IRAS saw Ganymed only twice and in both 
cases only the 25 /im flux was useable for the radiomet- 
ric calculations, resulting in a diameter of 31.7 ± 2.8 km 
and an albedo of 0.29 ± 0.06. A single chord occultation 
measurement gave a lower diameter limit of 16 km (Lan- 
gans 19850). The large lightcurve amplitude of 0.45 mag 



(Lagerkvist et al. |1989[ ) adds more uncertainties to the 



model calculations. Purely based on the given diameter 
and albedo, the TPM predicts approximately 0.1 Jy for 
the time of the ISOSS observation, which is well below 
the detection limit of this observing mode. 

5.2.3. No IRAS Detection 

5 asteroids (6 hits) have no IRAS detection, but fulfilled 
the conservative flux requirements for the ISOSS asteroid 
search (see Table ||). Unfortunately 4 sources (Euterpe, In- 
geborg, Cruithne and 1991 CS) have only marginal ISOSS 
detections and establishing upper flux limits was difficult. 
The results of Table | per Object: 



Metis: Kristcnscn ( 1984 ) has determined a size of 190 
± 19 km for this asteroid from an occultation event. A 
second occultation a few years later gave a high quality 
173.5 km diameter (Stamm |1989| , Blow |l997| ). Recent HST 
images (Storrs et al. 1999) revealed an elongated disk with 
a long axis of 235 km and a short axis of 165 km, which 
corresponds to an effective diameter of 197 km. Given 
the uncertainties involved we adopt the occultation result 
which is perfectly consistent with both techniques (see also 
Lagerros et al. 1999). The full light curve and shape in- 
formation has been taken from Erikson ( p000| ). The TPM 
predictions gave 4.1 ± 0.8 Jy and 3.4 ± 0.7 Jy, respec- 
tively (see Table ^|). Adopting the HST results instead led 
to about 5 % and 10 % higher fluxes. 

Based on the ISOSS flux of 3.8 Jy, the TPM allowed 
the calculation of an effective projected diameter of 178 km 
and an albedo of pv — 0.15 at the epoch of the ISOSS 
observation. A possible 20 % ISOSS flux uncertainty would 
correspond to about 10 % diameter uncertainty, resulting 



in a size of the rotating ellipsoid of 213 x 164 x 132 km 
with 10 % minimum uncertainties. 

Within the different uncertainties and based on the 
shape and spin vector solutions, the results agree nicely. 
The 3 methods -occultation, HST direct imaging and 
ISOSS radiometric method- led to comparable diameters 
and albedos. 



Euterpe: No IRAS observations are available. We used 
instead the largest extension from an occultation mea- 



surement (Dunham 1998) together with a shape and spin- 
vector model (Erikson |2000 ), H, G values (Piironen et al. 



1997 ) and an albedo of 0.13 related to the occultation 



cross section. The TPM prediction was 1.5 Jy with a large 
uncertainty due to the limited size knowledge. This is well 
within the detection limits, but the source was too far 
from the slew center to determine an upper flux limit. 

Ingeborg: There exists hardly any information about this 
asteroid. Based on its H-value of 10.1 mag, together with 
a typical S-type (Tholcn |1989| ) albedo of 0.155 one can 
calculate an approximate diameter of 32 km. The corre- 
sponding flux calculation for the ISOSS epoch gave 0.2 Jy, 
which is clearly below the detection limit. Even under the 
assumption of an extreme albedo of 0.03 the asteroid flux 
at 170 /jm would only be 1.3 Jy and therefore hardly de- 
tectable. Like in the case of Euterpe, Ingeborg had a slew 
center offset which was close to the maximal allowed 5'. 



Cruithne: Cruithne is currently the only known object 
on a horseshoe orbit around Earth (Christou 2000] ). It was 
also part of a special near-Earth object observations pro- 



gramme (Erikson et al. 2000a). Based on an unweighted 
mean of typical C and S-type asteroids (pv = 0.12) and 
an H-value of H = 15.13 ± 0.05, they calculated a diam- 
eter of 3.7 km and a slow rotation period of 27.4 hours. 
Although the observing geometry with only 0.37 AU from 
Earth was almost ideal, the 170 /im flux was only 0.1 Jy. 
Even an extremely low albedo (leading to a diameter of 
about 8 km) would only give 0.3 Jy well below the detec- 
tion limit. Therefore an upper limit from a background 
analysis would not give any new information. 



1991 CS: The case of 1991 CS is similar to Cruithne: 
A near-Earth asteroid, at only 0.14AU from Earth at 
the time of the ISOSS slew and with an H-value of 
17.4 mag. A radar campagne resulted in an estimated di- 
ameter of 1.1km, an albedo of 0.14 and a rotation pe- 



riod of 2.39 hours (Pravec et al. |1998| ). The TPM predicts 
a 170 /im flux below 0.1 Jy and even for extreme albedo 
values the flux would be below 0.3 Jy and therefore not 
detectable for ISOSS. 



|http : // sorry . vse . cz/~ludek/mp/world/mpoccl . txt 
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TDT 


Date/Time 


SSO 




-FModel 


Method 


Slew 


Real source 


Remarks 


JNo. 






(Jy) 


(Jy) 




speed 


offset* 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


83081700 


22-FEB-98 19:07:18 


(9) Metis 


3.8 


4.1 


1 


fast 


0' 


ok 


84281000 


06-MAR-98 16:31:55 


(9) Metis 


>1 


3.4 


2 


fast 


0.3' 


~ 3a detection 


63681300 


13-AUG-97 10:42:38 


(27) Euterpe 




1.5 




fast 


1.9' 


no detection 


25880900 


01-AUG-96 07:42:48 


(391) Ingeborg 




0.2 




fast 


2.5' 


no detection 


71682300 


01-NOV-97 06:03:56 


(3753) Cruithne 




0.1 




fast 


1.0' 


no detection 


27482100 


17-AUG-96 01:22:01 


(7822) 1991 CS 




0.1 




fast 


2.5' 


no detection 



Table 5. ISOSS observational results for asteroids without IRAS detection. * Closest approach to edge of closest pixel. 



5.2.4. Additional Results 

The Juno observations in Table |l| and in Table ^ have 
flux ratios systematically higher than ratios from compa- 
rable sources. Calibrating the ISOSS values with the cor- 
responding newly established methods 1 and 3a resulted 
in an average observation over model ratio of 1.14, in- 
dicating that the model diameter of Juno is about 7% 
too low. Miiller & Lagerros ( 2002a ) analysed 11 indepen- 
dent ISO observations, taken with the long wavelengths 
ISOPHOT detectors. They find a similar mean ratio of 
1.13±0.10, which confirms the tendency to higher diame- 
ter values. Both investigations indicate that the effective 
diameter should be close to 260 km, compared t o the pub- 
lished values of 241.4 km (Mii ller fc Lagerros |199S| ) and 
233.9±11.2km (Tedesco et al. [1992[ ). 

For Vesta the situation is not that clear. The values in 
Table |l| are not conclusive since the Vesta fluxes cover the 
difficult transition region between little flux loss and the 
more than 40% flux loss for sources brighter than 25 Jy 
(see Fig. |j). It seems that two of the measured fluxes 
(TDT 07881200, 79781500) are higher than one would 
expect from other sources with similar brightness. This 
contradicts the findings by Redman et al. ( |1998| ; |1992| ). 
They state for Vesta an extremely low emissivity of 0.6 in 
the submillimetre. Assuming that the emissivity is already 
lower in the far-IR would mean that the Vesta points in 
Fig. [| should lie below the general trend and not above. 
The measurement from methods 3a (TDT 57581500) and 
3b (TDT 61580800) agree within the errorba rs with the 
model predictions. As in Miiller & Lagerros ( E002aD , we 
see no clear indications of far-IR emissivities lower than 



the default values given in Miiller & Lagerros (1998) 



5.3. Comets 

5.3.1. Observational Results 

The results of the positional search through the ISOSS 
pointing data, combined with the flux estimates are given 
in Table The table columns are: (1-3) same as in Ta- 
ble @, (4-5) ISO-centric coordinates (2000.0), (6-7) Sun 



and Earth distance at the time of the observation, (8-9) 
ISOSS and model flux. 



2P (Encke): The comet has been detected at a slew 
end position on an extremely high background close to 
the galactic plane. The flux increase towards the comet 
nucleus corresponds to about 5-10 Jy. The coma exten- 
sion and its brightness profile could not be determined 
due to the high background brightness. The model flux at 
this close encounter with Earth (0.26 AU) may have been 
strongly overestimated due to the large apparent size of 
the central coma which was assumed to be of constant 
brightness. 

22P (Kopff): The ISOSS slew ends again on the comet, 
but this time the source is located on a clean low back- 
ground. The signal pattern is similar to that of a point- 
source with ~ 0.5-1 Jy, which is close to the detection 
limit. An upper flux limit of 2 Jy can be given, which is in 
good agreement with the simple model calculations (Ta- 
ble |). 

96P (Machholz 1): The comet has not been detected. 
The position calculation showed that the source was just 
outside the slew path, but within the specified 5' search 
limit. The low model flux indicated already the difficulty 
to detect the coma or the comet nucleus. 



103P (Hartley 2): Only a poor detection of an extended 
source was found, although the comet was on a low back- 
ground. ISOPHOT observations close to the ISOSS ob- 
serving epoch show that Hartley 2 had a colour tempera- 
ture of 285 K (Colangeli et al. |1999| ). This is 30 K colder 
than the calculated model temperature. A second reason 
for the discrepancy between a low ISOSS flux and the pre- 
dicted 25 Jy is probably the too large apparent size of the 
central coma which was assumed to be of constant bright- 
ness. At an Earth distance of 0.82 AU the model comet 
core covers a significant part of the aperture. Both effects 
together might explain the model value. 
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TDT 


Date/Time 
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RA. 






Dec. 






r 


A 


^Obs 


-Fmo 


IN (J. 








(hms^ 






(dms) 




( ATT"! 


l AU J 


\ J y) 


(•>y) 


(1) 




(2) 


(3) 


(4) 








(5) 






(6) 


(7) 


(s) 


(9) 


60780100 


14-JUL-97 22:57 


49 


2P/Encke 


14 


56 


28 


6 


-63 


36 


09 


1.164 


0.264 


5-10 


>50 


34881300 


29-OCT-96 22:27 


01 


22P/Kopff 


21 


23 


56 


6 


-19 


53 


22 


1.961 


1.523 


0.5-1 


<1 


23380800 


07-JUL-96 11:31 


51 


96P/Machholz 1 


23 


10 


27 


5 


-68 


19 


46 


2.052 


1.328 


no det. 


«1 


( ( ( oUZUU 


OL-UHi\j-SJ I 1U.OO 


DO 


1 fl^P /ttarf lo^ 9 

luor / ndi Liey z 




07 




1 

1 


07 




00 




89^ 


poor (let. 


~ZO 


80280100 


25-JAN-98 15:07 


48 


104P/Kowal 2 


Of) 


43 


46 





+08 


38 


35 


1.451 


1.496 


1 


w2 


33280100 


13-OCT-96 18:43 


22 


126P/IRAS 


21 


38 


46 


7 


-29 


54 


48 


1.712 


1.028 


— 


>2 


36280400 


12-NOV-96 13:51 


16 


126P/IRAS 


21 


45 


50 


3 


-08 


47 


05 


1.709 


1.307 


1.0 


<2 


13481800 


30-MAR-96 17:26 


33 


C/1995 Ol = 


19 


42 


20 


5 


-19 


43 


10 


4.867 


5.004 


ggg t6xt 




16280600 


27-APR-96 14:17 


19 


Hale-Bopp 


19 


44 


35 


1 


-17 


37 


42 


4.588 


4.259 


9.3+1.8 




31580500 


27-SEP-96 0:05 


16 




17 


29 


43 





-05 


11 


32 


2.934 


2.965 


30.9+7.3 




32081300 


01-OCT-96 17:21 


44 




17 


29 


42 


8 


-04 


57 


58 


2.878 


2.987 


see text 




32280200 


03-OCT-96 15:16 


59 




17 


29 


50 


7 


-04 


52 


28 


2.856 


2.995 


see text 




32580600 


06-OCT-96 23:28 


20 




17 


30 


13 


8 


-04 


42 


47 


2.816 


3.009 


see text 




77081500 


25-DEC-97 0:18 


51 




06 


32 


55 


7 


-64 


09 


08 


3.851 


3.683 


43.8+4.0 




86880300 


01-APR-98 14:02 


15 




05 


02 


13 


8 


-53 


09 


12 


4.855 


4.945 


see text 




87380400 


06-APR-98 16:12 


29 




05 


05 


07 





-52 


34 


21 


4.905 


5.009 


15.0+2.9 





Table 6. Observational geometry for the comets. No model values have been calculated for C/1995 Ol (Hale-Bopp). 
All hits are discussed in the text. One predicted hit was a "no detection", one a "poor detection" and in one case ( — ) 
the slew length was too short. 



104P (Kowal 2): A source of approximately 1 Jy was de- 
tected by one pixel at the predicted position of the comet, 
but confusion with a close IRAS source could not be ex- 
cluded. 

126P (IRAS): In the first case the slew length was only 1' 
which was not sufficient for the data analysis. The second 
case was a clear detection by one pixel (Method 2). The 
derived flux of 1 Jy is in agreement with the calculated 
upper limit of 2 Jy. 

C/1995 Ol (Hale-Bopp): 13481800: The slew passed 
over the coma with the nucleus only 30" outside the clos- 
est pixel. A weak signal of ^2MJysr _1 has been de- 
tected in this pixel. 16280600: The slew ended on Hale- 
Bopp and the integrated 4-pixel flux was determined to 
9.3±1.8Jy (Method 3a). 31580500: Method 3a was ap- 
plicable and a 170 /im flux of 30.9±7.3Jy has been de- 
rived. 32081300: Slew over the comet nucleus, with one 
pixel crossing centrally, two pixels in 1' distance and one 
pixel in 2' distance. The slew crossed the nucleus under 
an angle of 45° relative to the orientation of the dust 
tail (PsAngf\= 87.8°, ISOSSPosAng^ 42.0°). The mea- 



4 PsAng: The position angle of the extended Sun — » target 
radius vector as seen in the observer's plane-of-sky, measured 
counter-clock wise from reference frame North Celestial Pole. 

5 ISOSSPosAng: The position angle of the ISOSS slew 
origination as seen in the observer's plane-of-sky, measured 
counter-clock wise from reference frame North Celestial Pole. 



sured brightness profile clearly deviates from that of a 
point-source. The asymmetric profile is stronger towards 
the east, i.e. on the tail-side of the nucleus. Due to a nearby 
cirrus ridge, the dust tail extension is confirmed out to 2' 
only (but would be probably larger on a flat background) . 
32280200: The slew crossed the dust tail of Hale-Bopp 
under an angle of approx. 30° in 8' distance from the nu- 
cleus (PsAng = 87.1°, ISOSSPosAng = 55.8°). A signal 
increase at the position of the dust tail can be seen, but 
the signal pattern is difficult to discriminate from the cir- 
rus structures in the background, hence quantitatively not 
helpful. At the closest comet approach of 4' a second sig- 
nal increase can be seen which coincides with the position 
angle of the negative of the target's heliocentric velocity 
vector (Psj4MV[]= 151.2°). The signal increase is either 
related to the large cometary coma at a distance of only 
r = 2.86 AU from the sun or a kind of trail formation in the 
direction of PsAMV similar to what Reach et al. ( 2000 ) 
found for comet Encke. 32580600: The detectors moved 
centrally along the dust tail and cross over the comet nu- 
cleus (PsAng = 85.9°, ISOSSPosAng = 85.0°). The 
measured brightness profile clearly deviates from a point- 
source profile (see Figure ^) . A dust tail extension of more 
than 4' can be seen where the satellite approaches the 
nucleus. The signal in anti-tail direction decreases more 



rapidly (see also Section p. 3. 2j ). 77081500: Method 3a was 
applicable again and a 170 /im flux of 43.8±4.0 Jy was de- 



6 PsAMV: The position angle of the extended Sun — > target 
radius vector as seen in the observer's plane-of-sky, measured 
counter-clock wise from reference frame North Celestial Pole. 
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rived. 86880300: The slew passes in 4' distance ahead of 
the comet tail under an angle of approx. 70° with the sun 
direction (PsAng = 113.3°, ISOSSPosAng = 183.5°). A 
signal change of 2MJysr _1 extended over 15' can clearly 
be seen. Due to the viewing geometry (the phase an- 
gle is only 11°) coma and tail are difficult to separate 
and the signal increase is most likely connected to the 
dust emission of the extended coma and tail structures 
of Hale-Bopp. Here, as in 32280200, the signal maximum 
coincides with the Pa AMY angle of 1.2°. A connection 
might be possible between the 170 /^m signal pattern and 
large particles forming an elongated structure behind the 
comet nucleus while it is moving away from perihelion. 
87380400: Method 3a was applicable again and a 170 /im 
flux of 15.0±2.9 Jy was derived. 

5.3.2. C/ 1995 01 (Hale-Bopp) 



> 
1 
Pi 




-5 5 

Offset (arcmin) 

Fig. 6. An ISOSS profile of Hale-Bopp in comparison with 
the signal pattern of a point source. The slewing speeds 
for these two hits were identical and already quite low 
(2'/s), resulting in a high sampling rate. The offsets are 
defined with respect to the positions from the ephemeris 
calculations and have not been shifted for the two sources 
relative to each other. Hale-Bopp shows an asymmetri- 
cally extended profile. Note that the profiles are not de- 
convolved. 

Figure ^| illustrates a measured signal profile from a 
central slew over Hale-Bopp (TDT 32580600) in compar- 
ison with a slew over the point source Neptune (TDT 
72081600). For these two detections, the geometrical con- 
figurations and the slewing speeds (2'/s) have been identi- 
cal. The detectors moved first centrally over the dust tail 
(left side of the peak) and then over the nucleus of Hale- 
Bopp (peak). An asymmetric signal profile between —6' 
and +4' can be seen. In case of Neptune the signal in- 
crease starts approximately 2' ahead of the true position, 



which is related to a combination of the Airy disk with 
the slew speed and read out frequency (see also Hotzel et 
al. |200l| ). The slight shift between the two peaks is most 
probably related to the possible positional uncertainties of 



ISOSS data (see Section 4.4). This Hale-Bopp asymmetry 
has not been seen in dedicated 170 /im maps (Peschke et 
al. 1999) which were taken at r = 3.904 AU (as compared 
to r = 2.816 AU in Figure |). 

The fluxes derived fro m Me thod 3a can be compared 
to results from Grim et al. 2001 through the following cor- 
rections: 1) Flux corrections according to Figure^; 2) Nor- 
malization to a standard aperture diameter of 23" assum- 
ing that the coma brightness scales linearly with aperture 
diameter (c a 23 = 0.120); 3) Point-spread- function correc- 
tion which takes into account the differences of a point 
source and a l/p-coma (c ps f — 1.092); 4) Colour correc- 
tion^ which changes for different distances from the sun 

{^colour — fv))' 

The first two measurements (16280600, 31580500 ) were 
obtained on the same days as the ones in Grim et al. |2001| 



The calibrated and reduced 23" fluxes agree within the 
errorbars. The third observation (77081500), taken 5 days 
earlier than the dedicated Hale-Bopp observation, lead to 
a flux of 5.64±0.56Jy, as compared to 2.66 Jy. This large 
difference can not be explained by epoch or geometry dif- 
ferences. However, the dedicated measurement was mis- 
pointed by 24" which caused large uncertainties in the ap- 
plied corrections. The ISOSS flux provides therefore valu- 
able information for the colour temperature determination 
and, through grain size models, might give clues whether 
icy grains were present in the coma in December 1997 
at almost 4AU post-perihelion. The last measurement of 
Method 3a (87380400) was obtained when Hale-Bopp was 
already 4.9 AU from the sun. The calibrated and reduced 
23" flux was 1.97±0.39 Jy. This is the most distant thermal 
far-IR observation of Hale-Bopp post-perihelion. A com- 
parison of the flux with a dedicated observation (Grim 
ct al. 



2001 



F v = 1.06 Jy) at a similar distance from the 
sun, but pre-perihelion, shows that the dust emission post- 
perihelion was higher as the comet receeded from the sun. 
In fact, the higher far-IR fluxes post-perihelion are related 
to contributions from large particles which have been ac- 
cumulated during the passage around the sun and which 
stay on similar orbits as the nucleus. 

Two measurements (32280200, 86880300) show signal 
peaks a few arcminutes away from the nucleus in anti- 
orbital velocity (trail) direction. It seems that the emit- 
ting dust particles are not homogeneously distributed and 
are concentrated in a narrow region of the outer parts of 
the dust coma towards the trail direction. These features 
are not seen in slews over other parts of the outer coma. 



7 There seems to be a wrong application (multiplication in- 
stead of division) of the colour correction factor by Grim et al. 



(2001). For consistency, we do however apply all corrections as 



given in their paper. 
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Reach et al. 2000 observed for the first time the dust trail 



formation in comet Encke in the mid-IR. The ISOSS data 
provide now evidence for this process in the far-IR where 
the emission is strongly connected to the largest particles. 



6. Conclusions and Outlook 



The purpose of the SSO exlractioii from the ISOSS was 

manifold: Calibration aspects, catalogue cleaning aspects 
and scientific aspects. The main achievements were clearly 
in the calibration section, were serendipitously seen as- 
teroids and planets led to an improved flux calibration 
for ISOSS targets. Bright sources of the automatic point- 
source extraction procedure have now a solid calibration 
basis. It was also possible to establish new methods to 
calibrate source detections under a large variety of cir- 
cumstances, including the important slew end positions 
and low slewing speeds. 

The aspect of SSO cleaning from ISOSS catalogue lists 
will avoid wrong identifications and help follow up pro- 
grammes of galactic and extra-galactic sources. 

The outcome of the scientific analysis of SSO detec- 
tions were modest due to the limitations of the ISOSS 
mentioned in Sect. [| Despite all difficulties we could 
demonstrate that the far-IR fluxes of asteroids are impor- 
tant. Diameter and albedo estimates through TPM cal- 
culations are much more reliable than estimates based on 
visible brightness alone. An accurate H- value of 12.0 mag 
would allow for diameters ranging from 10.6 km (pv = 
0.25) to 23.7km (p v = 0.05), corresponding to a ±40% 
uncertainty of the average. An additional thermal flux 
with a flux error of ±20 % allows a 4 times more accu- 
rate diameter determination. 

The ISOSS results for Hale-Bopp are more valuable. 
They can now be used for additional comet modeling (e.g. 
models by Hanner 1983) for more reliable interpretation of 
grain properties and ice influences at different heliocentric 
distances. Comets are usually optically bright due to fresh 
ice surfaces, but in the far-IR the sublimated larger par- 
ticles dominate the thermal emission. After many orbits 
around the sun these large particles form trails which were 



first measured by IRAS (Sykes 1986). For Hale-Bopp we 



found significantly more thermal emission post-perihelion 
than for comparable configurations pre-perihelion. Addi- 
tionally we saw asymmetries due to the dust tail and an in- 
dicative detection of large particles concentrated towards 
the anti-orbital velocity (trail) direction. 

The expectations for future far-IR and submillimetre 
projects on SSO related topics are large: SIRTF, SOFIA, 
ASTRO-F, HERSCHEL and others will have many dedi- 
cated programmes on asteroids, comets and planets, but 
will also see by chance interesting targets. Especially the 
ASTRO-F/FIS all sky survey in 4 photometric bands in 
the region 50 to 200 /im will serendipitously detect many 
SSOs. Our experience with ISOSS in terms of calibration 



through asteroids and planets, but also in identification of 
moving targets could then be of great benefit. 
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